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Abstract 



Preference handling and optimization are indispensable means for addressing non-trivial applications 
in Answer Set Programming (ASP). However, their implementation becomes difficult whenever they 
^^ . bring about a significant increase in computational complexity. As a consequence, existing ASP 

systems do not offer complex optimization capacities, supporting, for instance, inclusion-based min- 
imization or Pareto efficiency. Rather, such complex criteria are typically addressed by resorting to 



"^ ■ dedicated modeling techniques, like saturation. Unlike the ease of common ASP modeling, however, 

' these techniques are rather involved and hardly usable by ASP laymen. We address this problem by 

developing a general implementation technique by means of meta-programming, thus reusing exist- 
ing ASP systems to capture various forms of qualitative preferences among answer sets. In this way, 

^ ' complex preferences and optimization capacities become readily available for ASP applications. 

lO ' 1 Introduction 



Preferences are often an indispensable means in modeling since they allow for identi- 
fying preferred solutions among all feasible ones. Accordingly, many forms of prefer- 
ences have already found their way into systems for Answer Set Programming (ASP; 
(IBaral 20031) ). For instance, smodels provides optimization statements for expressing cost 
functions on sets of weighted literals ( iSimons et al. 2002] ). and dlv (Leone e t al. 2006l l of- 
fers weak constraints for the same purpose. Further approaches ( ,Delgrande et al. 2003 [ 
C^ \ lEiter et al. 20031) allow for expressing various types of preferences among rules. Unlike 

this, no readily applicable implementation techniques are available for qualitative pref- 
erences among answer sets, like inclusion minimahty, Pareto-based preferences as used 
in (ISakama and Inoue 20001 IBrewka et al. 2004l l. or more complex combinations as pro- 
posed in (IBrewka 20041) . This shortcoming is due to their higher expressiveness leading to 
a significant increase in computational complexity, lifting decision problems (for normal 
logic programs) from the first to the second level of the polynomial time hierarchy (cf. 



( Garey and Johnson 1979 1). Roughly speaking, preferences among answer sets combine 
an NP with a coNP problem. The first one defines feasible solutions, while the second 
one ensures that there are no better solutions according to the preferences at hand. For 
implementing such problems, Eiter and Gottlob invented in (119951 1 the saturation tech- 
nique, using the elevated complexity of disjunctive logic programming. In stark contrast 
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to the ease of common ASP modeling (e.g., strategic companies can be "naturally" en- 
coded jLeone et al. 2006 ) in disjunctive ASP), however, the saturation technique is rather 
involved and hardly usable by ASP laymen. 

For taking this burden of intricate modeling off the user, we propose a general, 
saturation-based implementation technique capturing various forms of qualitative prefer- 
ences among answer sets. This is driven by the desire to guarantee immediate availabiUty 
and thus to stay within the realm of ASP rather than to build separate (imperative) com- 
ponents. To this end, we take advantage of recent advances in ASP grounding technology, 
admitting an easy use of meta-modeling techniques. The idea is to reinterpret existing opti- 
mization statements in order to express complex preferences among answer sets. While, for 
instance in smodels, the meaning of ^minimize is to compute answer sets incurring min- 
imum costs, we may alternatively use it for selecting inclusion-minimal ones. In contrast 
to the identification of minimal models, investigated by Janhunen and Oikarinen in (120041 
I2OO8I 1. a major challenge lies in guaranteeing the stability property of implicit counterex- 
amples, which must be more preferred answer sets rather than (arbitrary) models. For this 
purpose, we develop a refined meta-program qualifying answer sets as viable counterex- 
amples. Unlike the approach of Eiter and Polleres (120061 ). our encoding avoids "guessing" 
a level mapping to describe the formation of a counterexample, but directly denies mod- 
els for which there is no such construction. Notably, our meta-programs apply to (reified) 
extended logic programs (ISimons et al. 2002] l. possibly including choice rules and #suni 
constraints, and we are unaware of any existing meta-encoding of their answer sets, neither 
as candidates nor as counterexamples refuting optimality. 



2 Background 

We consider extended logic programs (ISimons et al. 2002] l allowing for (proper) disjunc- 
tions in heads of rules dGelfond and Lifschitz 199fl l. A rule r is of the following form: 

H <— Bi, . . . , Bm, ^Bm+l, • • ■ , ^Bn- 

By head{r) = H and hody{r) = {Si, . . . , i?m, ^Bm+i, ■ ■ ■ , ^Bn}, we denote the head 
and the body of r, respectively, where "~" stands for default negation. The head H is a 
disjunction ai V • • • V a/c over atoms ai, . . . ,ak, belonging to some alphabet A, or a ^sum 
constraint L ^sum.[ii = wi, . . . ,ik = Wk] U. In the latter, £i = Ui or £i — ^ai is a literal 
and Wi a non-negative integer weight for ai ^ A and 1 < i < k; L and U are integers 
providing a lower and an upper bound. Either or both of L and U can be omitted, in which 
case they are identified with the (trivial) bounds and cxd, respectively. A rule r such that 
head{r) = _L (H is the empty disjunction) is an integrity constraint. Each body component 
Bi is either an atom or a ^^^sum constraint for 1 < i < n. If hody{r) = 0, r is called a. fact, 
and we skip"<— " when writing facts below. For aset {_Bi, . . . , Brm ^B„i+i, . . . , ^Bn}, a 
disjunction ai V • • ■ V a^, and a T^sum constraint L 7^sum[£i = wi, . . . ,£k = Wk] U, we let 
{Bi, . . . , Bm, ^Bm+i, ■ . ■ , ^Bn}^ == {Bi,...,Bm}, (ai V ••• Vafc)+ = {ai,...,afc}, 
and (L #suni[^i = wi, . . . ^Ik = Wk] U)+ = [ii ^ Wi\l <i <k,£^ <^ A\. Note that the 
elements of a #sum constraint form a multiset, possibly containing duplicates. For some 
S = {ai, . . . , Ofc} or S = [ai = wi, . . . , a^ = Wk], we define atom{S) = {ai, . . . , Uk}- 
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A (Herbrand) interpretation is represented by the set X C ^ of its entailed atoms. The 
satisfaction relation "|=" on rules r is inductively defined as follows: 

• X \^^BifX ^ B, 

• X h (ai V • • • V flfe) if {ai, . . . , afc} n X 7^ 0, 

• X \= (L#sum[£i =wi,...Jk^ Wk]U) if L < J2i<i<k.x^i,^^ ^ ^^ 

• X \= body{r) if X h £ for all £ G body{r), and 

• X \=rifX \= head{r) or X ^ hody{r). 

A logic program 11 is a set of rules r, and X is a model of 11 if X |= r for every r G 11. 
The reduct of the head H of a rule r wrt X is H^ ~ {ai V • • • Vofc} if iJ = ai V • • • Va^, 
and H^ = atom{H^) D X if H — L #sum[^i — wi, . . . ,ik — Wk] U . Furthermore, the 
reduct of some (positive) body element B e body{r)^ is B^ — B if B ^ A, and B^ = 
{L~J2i<i<kA=~ai.,a,<^x'^i)#^^'^B+ifB = L#sum[^i ^wi,...,£k^ Wk]U. The 
reduct of 11 wrt X is the following logic program: 

n'' = {H^B^,...,B^,\ 

r en,X 1= hody(r),H E head (r)^ , body {r)+ = {Bi, . . . , B„^}} . 

That is, for all rules r G 11 whose bodies are satisfied wrt X, the reduct is obtained by 
replacing #sum constraints in heads with individual atoms belonging to X and by elimi- 
nating negative components in bodies, where lower bounds of residual #sum constraints 
(with trivial upper bounds) are reduced accordingly. Finally, X is an answer set of 11 if X 
is a model of 11 such that no proper subset of X is a model of 11^. In view of the latter 
condition, note that an answer set is a minimal model of its own reduct. 

The definition of answer sets provided above applies to logic programs containing ex- 
tended constructs (#suin constraints) under "choice semantics" JSimons et al. 2002] l. while 
additionally allowing for disjunctions under minimal-model semantics (wrt a reduct). We 
use these features to embed extended constructs of an object program into a disjunctive 
meta-program, so that their combination yields optimal answer sets of the object program. 
To this end, we reinterpret =i^minimize statements of the following form: 

#minimize[^i = u;i@Ji, . . . , 4 = Wk@Jk]- (1) 

Like with #sum constraints, every £i is a literal and every Wi an integer weight for 
1 < i < fc, while Ji additionally provides an integer priority level]^ Priorities allow 
for representing a sequence of lexicographically ordered ^minimize objectives, where 
greater levels are more significant than smaller ones. By default, a ^minimize statement 
distinguishes optimal answer sets of a program 11 in the following way. For any X C A 
and integer J, let S^ denote the sum of weights w over all occurrences of weighted literals 
£ = w@J in ^ such that X ^ £. An answer set X of 11 is dominated if there is an answer 
set F of n such that Ej < E j and E j, = J2j, for all J' > J, and optimal otherwise. 
In the following, we assume that every logic program is accompanied with one (possibly 

^ Explicit priority levels are supported in recent versions of the grounder gringo jGebser et al.t . This avoids 
a dependency of priorities on input order, which is considered by Iparse I Syrjaneni if several ^minimize 
statements are provided. Priority levels are also supported by dlr (Leone et al. 2006i in weak constraints. Fur- 
thermore, we admit negative weights in ^minimize statements, where they cannot raise semantic problems 
(cf. jFerraris 2005) ) going along with the rewriting of #;sujii constraints suggested in l iSimons et al. 2002; . 
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empty) ^^minimize statement of the form ([T]i- Instead of the default semantics, we con- 
sider Pareto efficiency wrt priority levels J, weights w, and several distinct optimization 
criteria. In view of this, we use levels for inducing a lexicographic order, while weights are 
used for grouping literals (rather than summation). Pareto improvement then builds upon 
a two-dimensional structure of orderings among answer sets, induced by J and w. In turn, 
each such pairing is associated with some of the following orderings. By Y < j X, we 
denote that the cardinality of the multiset of occurrences of I = w@J in ^ such that 

Y \= i is not greater than the one of the corresponding multiset for X \= £. Furthermore, 
we write Y C^ X if, for any weighted literal £ = ■w@J occurring in (HJ, Y \= ( implies 
X \^ £. As detailed in the extended version of this paper JGebser et al. 201 ll) . we addi- 
tionally consider the approach of JSakama and Inoue 2000l l and denote by Y ^ j X that 

Y is preferable to X according to a (given) preference relation ^ among literals £ such 
that £ = w@J occurs in ([U. Given a logic program 11 and a collection M of relations of 
the form Oj for priority levels J, weights w, and o G {<, C, ^}, an answer set y of 11 
dominates an answer set X of 11 wrt M if there are a priority level J and a weight w such 
that X o}" Y does not hold for o^ G M, while Y oj X holds for all oj G M where 
J' > J. In turn, an answer set X of 11 is optimal wrt M if there is no answer set y of 11 
that dominates X wrt M. 

As an example, consider the following program, referred to by Ho: 

(2) 
(3) 
(4) 

This program has five answer sets, viz. {p,q}, {p,r}, {p, s}, {p, s,t}, and {s,t}. (Sets 
{ai, . . . , flfc} in dU and (|3]l are used as shorthands for =i^sum[ai = 1, . . . , a^ — 1].) In 
addition, let Hi denote the union of Ho with the following ^minimize statement: 

#iniiiimize[p = l@l,q= l@l,r = l@l,s = 1@1]. (5) 

This statement specifies that all atoms of IIq except for t are subject to minimization. Pass- 
ing Hi to gringo and an answer set solver like smodels yields the single < \ -minimal answer 
set {s, i). Note, however, that IIq has three cj-minimal answer sets, namely {p, q], {p, r}, 
and {s, i). They cannot be computed directly from Hi via any available ASP system. 

We implement the complex optimization criteria described above by meta-interpretation 
in disjunctive ASP. For transparency, we provide meta-programs as true ASP code in the 
first-order input language of gringo (IGebser et al.l l. including not and | as tokens for ^ 
and V, respectively, as well as {ai , . . . , ak} as shorthand for #suin [ai = l , . . . , afc=l ] . 
Further constructs are informally introduced by need in the remainder of this paper Note 
that our (disjunctive) meta-programs apply to an extended object program that does not 
include proper disjunctions (over more than one atom). Unless stated otherwise, we below 
use the term extended program to refer to a logic program without proper disjunctions. 

3 Basic Meta-Modeling 

For reinterpreting ^^minimize statements by means of ASP, we take advantage of recent 
advances in ASP grounding, admitting an easy use of meta-modeling techniques. To be 
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1 rule (pos (sum (l,0,2)),pos(conjunction{0)}). % 1 { p, t } :- 1 { r, s, not t } 2. 

2 wlist (0, O,pos (atom(p) ) , 1) . wlist (0, 1, pos (atom (t) ), 1) . 

3 set (O,pos (sum(l, 1,2))). 

4 wlist (1, 0, pos (atom(r) ), 1) . wlist (1, 1, pos (atom (s) ), 1) . wlist ( 1, 2 , neg (atom (t )), 1 ) . 

6 rule (pos (sum (0,2,l)),pos(conjunction(l))). % { q, r } 1 :- 1 { p, t }. 

7 wlist (2, O,pos (atom(q) ) , 1) . wlist (2, 1, pos (atom (r) ), 1) . 

8 set (l,pos (sum(l, 0, 2) ) ) . 

10 rule(pos( at om(s)),pos(con junction (2))). % s :- not q, not r. 

11 set (2, neg (atom(q) ) ) . set (2, neg (atom(r) ) ) . 

13 scc(0,pos (atom (p))). scc(0,pos (atom (r))). scc(0,pos (atom (t ) ) ) . 

14 scc(0,pos(conjunction(0))). scc(0,pos (sum (1,1,2))). 

15 scc(0,pos(conjunction(l))). scc(0,pos (sum (1,0,2))). 

17 minimize (1,3). % #minimize [p-l(3 1, q-l(3 1, r-l@l, s-l@l]. 

18 wlist (3, O,pos (atom(p) ) , 1) . wlist (3, 1, pos (atom (q) ), 1) . 

19 wlist (3,2,pos (atom(r) ) , 1) . wlist (3, 3, pos (atom (s) ) , 1) . 



Listing 1 . Facts describing a reified extended logic program. 



precise, we rely upon the unrestricted usage of function symbols and program reification 
as provided by gringo (iGebser et al.l) . The latter allows for turning an input program along 
with a ^minimize statement into facts representing the structure of their ground instanti- 
ation. 

For illustrating the format output by gringo, consider the facts in Line 1-15 of Listing[Tl 
obtained by calling gringo with option --reify on program Ho. Let us detail the repre- 
sentation of the rule in (|2|i inducing the facts in Line 1-4. The predicate rule/ 2 is used to 
link the rule head and body. By convention, both are positive rule elements, as indicated via 
the functor pos /I. Furthermore, the term sum (1, 0,2) tells us that the head is a T^sum 
constraint with lower bound 1 and (trivial) upper bound 2 over a list labeled of weighted 
literals. In fact, the included literals are provided via the facts over wlist/ 4 given in 
Line 2, whose first arguments are 0. While the second arguments, and 1, are simply 
indexes (enabling the representation of duplicates in multisets), the third ones provide lit- 
erals, p and t, each having the (default) weight 1, as given in the fourth arguments. Again 
by convention, the body of each rule is a conjunction, where the term conjunction (0) 
in Line 1 refers to the set labeled 0. Its single element, a positive #sum constraint with 
lower bound 1 and upper bound 2 over a list labeled 1, is provided by the fact in Line 3. 
The corresponding weighted literals are described by the facts in Line 4; observe that the 
negative literal not t is represented in terms of the functor neg/ 1, applied to atom (t ) . 
The rules in (O and (HI are represented analogously in Line 6-8 and 10-11, respectively. It 
is still interesting to note that recurrences of lists of weighted literals (and sets) can reuse 
labels introduced before, as done in Line 8 by referring to . In fact, gringo identifies repe- 
titions of structural entities and reuses labels. In addition to the rules of Hq, the elements of 
non-trivial strongly connected components of its positive dependency graph (cf . Q below) 
are provided in Line 13-15. Albeit their usage is explained in the next section, note already 
that the members of the only such component, labeled 0, include atoms as well as (posi- 
tive) body elements, i.e., conjunctions and #sum constraints, connecting the component. 



Indeed, the existence of facts over scc/2 tells us that Hq is not tight (cf. (Fages 1994l). 
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1 % extract rule elements 



3 lltb(B) 

4 lltb(E) 

5 lltb(E) 



rule (_,B) . 

litb (pos (conjunction (S) } ) , set {S, E) . 

eleb{sum(_, S,_) ) , wllst (S,_,E,_) . 



7 eleb(P) :- litb (pos (P) ) . 

8 eleb(N) :- litb (neg (N) ) . 



10 elem(E) 

11 elem(E) 

12 elem(P) 

13 elem(N) 



eleb(E) . 

rule (pos (E) ,_) . 

rule (pos (sum (_, S, _)),_) , wlist (S,_, pos (P) ,_) . 

rule (pos (sum (_, S, _)),_) , wlist (S,_,neg(N) ,_) . 



15 % generate answer set from reified rules 

17 hold (con junction (S) } :- eleb (con junction (S) ) , 

18 hold(P) : set (S,pos (P) ) , 

19 not hold(N) : set (S, neg (N) ) . 

20 hold(sum(L, S,U) ) :- eleb (sum (L, S, u) ) , 

21 L #sum [ hold(P) = W : wlist (S, Q, pos (P) , W) , 

22 not hold(N) = W : wlist (S, Q, neg (N) , W) ] U. 

24 hold(atom(A) ) :- rule (pos (atom (A) ) , pos (B) ) , hold(B). 

25 L #sum [ hold(P) = W : wlist (S, Q, pos (P) ,W) , 

26 not hold(N) = w : wlist (S, Q, neg (N) ,W) ] u 

27 :- rule(pos (sum(L, S,U) ) ,pos (B) ) , hold(B). 

28 :- rule (pos (false) , pos (B) ) , hold(B). 

30 % project output to atoms of answer set 
32 #hide. #show hold (atom (A) ) . 

Listing 2. Basic meta-program (meta . Ip) for reified extended logic programs. 

Now, we may compute all five answer sets of Hq (given in pO . Ip) by combining the 
facts in Line 1-15 of Listing [T] with the basic meta-program in Listing |2] (met a . lp)o 

gringo — reify pO.lp | gringo meta.lp - | clasp 

Each answer set of the meta-program applied to a reified program corresponds to an an- 
swer set of the reified program. More precisely, a set X of atoms is an answer set of 
the reified program iff the meta-program yields an answer set Y such that X ~ {a \ 
hold (atom (a) ) G F}, e.g., hold (atom (q) ) stands for q. As indicated in the com- 
ments (preceded by %), our meta-program consists of three parts. Among the rule elements 
extracted in Line 3-13, only those occurring within bodies, identified via eleb/1, are 
relevant to the generation of answer sets specified in Line 17-28. (Additional head ele- 
ments, given by elem/1, are of interest in the next section.) In fact, answer set generation 
follows the structure of reified programs, identifying conjunctions and ^^sum constraints 
that holqfl to further derive atoms occurring in rule heads, either singular or within ^/^sum 
constraints (cf. Line 24-27). Line 28 deals with integrity constraints represented via the 
constant false in heads of reified rules. The last part in Line 32 restricts the output of the 
meta-program's answer sets to the representations of original input atoms. 

Finally, note that meta . Ip does not inspect facts representing a reified ^minimize 
statement, such as the ones in Line 17-19 of Listing[T]stemming from the statement in (|5]l. 



^ Following Unix customs, tlie minus symbol "-" stands for tlie output of "gringo — reify pO . Ip." 
^ The " ; " connective expands to the Hst of all instances of its left-hand side such that corresponding instances 
of literals on the right-hand side hold (cf. (iSyrjanenj and liGebser et al.,!). 
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Such facts over minimize/2 provide a priority level as the first argument and the label of 
a list of weighted literals, like the ones referred to from within terms of functor sum/ 3, as 
the second argument. Rather than simply mirroring the standard meaning of ^minimize 
statements (by encoding them analogously to rules; cf. Line 17-28 of Listing|2l), we support 
flexible customizations. In fact, the next section presents our meta-programs implementing 
preference relations and Pareto efficiency, as described in the background. 

4 Advanced Meta-Modeling 

Given the reification of extended logic programs and the encoding of their answer sets in 
met a . Ip, our approach to complex optimization is based on the idea that an answer set 
generated via met a . Ip is optimal (and thus acceptable) only if it is not dominated by any 
other answer set. For implementing our approach, we exploit the capabilities of disjunctive 
ASP to compactly represent the space of all potential counterexamples, viz. answer sets 
dominating a candidate answer set at hand. To this end, we encode the subtasks of 

1 . guessing an answer set as a potential counterexample and 

2. verifying that the counterexample dominates a candidate answer set. 

A candidate answer set passes both phases if it turns out to be infeasible to guess a coun- 
terexample that dominates it. For expressing the non-existence of counterexamples, we 
make use of an error-indicating atom bot and saturation jEiter and Gottlob 19951) . deriv- 
ing all atoms representing the space of counterexamples from bot. Since the semantics 
of disjunctive ASP is based on minimization, saturation makes sure that bot is derived 
only if it is inevitable, i.e., if it is impossible to construct a counterexample. However, via 
an integrity constraint, we can stipulate bot (and thus the non-existence of counterexam- 
ples) to hold, yet without providing any derivation of bot. In view of such a constraint 
and saturation, a successful candidate answer set is accompanied by all atoms representing 
counterexamples. Given that the reduct drops negative literals, the necessity that all atoms 
representing counterexamples are true implies that we cannot use their default negation in 
any meaningful way. Hence, we below encode potential counterexamples, i.e., answer sets 
of extended programs, and (non-)dominance of a candidate answer set in disjunctive ASP 
without taking advantage of default negation (used in met a . Ip). 

For encoding the first subtask of guessing a counterexample, we rely on a characteriza- 
tion of answer sets in terms of an immediate consequence operator T (cf. (Lloyd 1987[ )), 
defined as follows for a logic program 11 and a set X C ^ of atoms: Tu{X) = {head{r) \ 
r G n, X ^ body{r)}. Furthermore, an iterative version of T can be defined in the fol- 
lowing way: T-^iX) = X and T^^^^iX) = Tt{{X) U rn(Ti^(^)). In the context of an 
extended program 11, possibly including choice rules, default negation, and upper bounds 
of weight constraints, we are interested in the least fixpoint of T applied wrt the reduct n"^ . 
Since a fixpoint is reached in at most | atom (11) | applications of T, where atom (11) C A 
denotes the set of atoms occurring in 11, the least fixpoint is given by Tjjx"™ (0)- As 
pointed out in dLiu and You 201 Oi l, a model X of an extended program 11 is an answer 
set of n iff rjix°"^"^'(0) = X. Furthermore, Liu and You (l20T0l l show that X violates 
the loop formula of some atom or loop if X is a model, but not an answer set of 11. This 
property motivates a "localization" of T on the basis of (circular) positive dependencies. 
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The (positive) dependency graph of an extended program 11 is given by the following 
pair of nodes and directed edges: 

(atom(n),{(a, &) | r £ n,a £ atom {head (r)^ ), B £ body{r)^,b £ atoTn{B^)}) . (6) 

A strongly connected component (SCC) is a maximal subgraph of the dependency graph 
of n such that all nodes are pairwisely connected via paths. An SCC is trivial if it does not 
contain any edge, and non-trivial otherwise. Note that the SCCs of the dependency graph 
of n induce a partition of atom (11) such that every atom and every loop of 11 is contained 
in some part. Hence, we can make use of the partition to apply T separately to each part. 

Proposition 1 

Let n be an extended logic program, Ci, . . . ,Ck be the sets of atoms belonging to the 

SCCs of the dependency graph of 11, and X C atom (11). 

Then, we have that Tjix""'"^' (0) = ^ iff Ui<j<fc(Tnx ' (^ \ Q) n Q) = X. 

The above property is used in our encoding of answer sets (as counterexamples) in dis- 
junctive ASP. In a nutshell, it combines the following parts: 

1 . guessing an interpretation, 

2. deriving the error-indicating atom bot if the interpretation is not a supported model 
(where each true atom occurs positively in the head of some rule whose body holds), 

3. deriving bot if the true atoms of some non-trivial SCC are not acyclicly derivable 
(checked via determining the complement of a fixpoint of T), and 

4. saturating interpretations that do not correspond to answer sets by deriving all truth 
assignments (for atoms) from bot. 

Note that the third part, checking acyclic derivability, concentrates on atoms of non-trivial 
SCCs, while checking support in the second part is already sufficient for trivial SCCs. 

The meta-program in Listing [3] implements the sketched idea. In the following, we con- 
centrate on describing its crucial features. For evaluating support, the meta-rules in Line 3 
and 4 collect atoms having a positive occurrence in the head of a rule along with the 
rule's body. Note that, for atoms contained in a #siim constraint in the head, the asso- 
ciated bounds and weights are inessential in the context of support. On the other hand, the 
meta-rule in Line 6 sums the weights of all literals in a #sum constraint; this is needed 
to evaluate bounds in the sequel, where (non-reified) default negation and upper bounds 
(acting negatively) are inapplicable in view of saturation. 

The meta-rules in Line 10-29 generate an interpretation by guessing some truth value 
for each atom (Line 10) and evaluating further constructs occurring in a reified program 
accordingly (Line 12-29). While the special constant false (used as head of integrity 
constraints) holds in no interpretation (fail (false) is a fact) and the evaluation of 
conjunctions is straightforward, more care is required for evaluating #sum constraints. For 
instance, the case that a ^^sum constraint holds is in the meta-rule in Line 19-23 identified 
via sufficiently many literals that hold to achieve the lower bound L and also sufficiently 
many hterals that do not hold to fill the gap between the upper bound U and the sum T of 
all weights. Note that the latter condition is encoded by the lower bound T-U, rather than 
taking U as an upper bound (as done in met a . Ip). The complementary cases that a #sum 
constraint does not hold are described in the same manner in Line 24-29, where the lower 
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% extract supports of atoms and sums of weight lists' weights 

supp (atom (A),B) :- rule(pos{atom(A)), pos{B)). 

supp (atom (A),B) :- rule(pos (sum (_,S,_)),pos(B)}, wlist(S,_, pos (atom (A) ) ,_) . 

sum(S,T) :- elem(sum(_, S,_) ) , T = #sum [ wlist (S, Q,_, W) = W ]. 

% generate interpretation 

true (atom (A) ) I fail (atom (A) ) :- elem (atom (A) ) . 

fail (false) . 



14 
15 
16 
17 

19 

20 
21 
22 
23 
24 
25 
26 
27 
28 
29 



true (conjunction (S) } 

fail (conjunction (S) } 
fail (conjunction (S) } 

true (sum(L, S,U) ) 



fail(sum(L,S,U) ) 



fail (sum(L, S,U) ) 



elem (conjunction (S) ) , 

true(P) : set (S,pos (P) ) , fail (N) : set (S, neg (N) ) . 
elem (conjunction(S}), set(S,pos(P)}, fail(P}. 
elem (con junction (S) ) , set (S, neg (N) } , true(N}. 



elem(sum(L, S,U) ) , sum(S,T), 



wlist (S,Q,pos (P) , W) , 
wlist (S,Q, neg (N) ,W) ], 
wlist (S,Q,pos (P) , W) , 
wlist (S,Q, neg (N) ,W) ]. 



L #sum [ true (P) = W 

fail (N) = W 

T-U #sum [ fail (P) = W 

true (N) = W 

:- elem(sum(L, S,U) ) , sum(S,T), 

T-L+1 #sum [ fail(P) = W : wlist (S, Q, pos (P) , W) , 

true(N) = W : wlist (S, Q, neg (N) , W) ]. 
:- elem(sum(L, S,U) ) , 

U+1 #sum [ true(P) = W : wlist (S, Q, pos (P) , W) , 

fail(N) = W : wlist (S,Q, neg (N) ,W) ]. 



31 % verify supported model properties 



33 
34 



36 

38 



40 

41 



bot :- rule (pos (H) , pos (B) ) , true (B) , fail(H). 
bot :- true (atom (A) ) , f ail (B) : supp (atom (A) , B) . 



% verify acyclic derivability 



step(C,Z) 
scow (A) 



43 wait(E,D-l) 



45 
46 

47 

49 
50 
51 
52 

54 
55 

57 

59 

61 
62 



wait (atom (A) , 0) 
wait (atom (A) ,D) 



scc(C,_}, Z - #sum [ sec (C, pos (atom (A) ) } ]. 

sec (C,pos (atom(A} ) } , 

fail (B) : supp (atom (A) , B) : not sec (C, pos (B) ) . 

sec (C,pos (E) ) , fail (E) , step(C,Z), D = 1..Z. 

sec (C, pos (atom (A) ) } . 

sec (C,pos (atom(A) ) ) , scow (A) , step(C,Z), D = 1..Z, 

wait(B,D-l) : supp (atom(A) ,B) : sec (C, pos (B) ) . 



wait (sum(L, S,U) ,D-1) :- sec (C, pos (sum (L, S, U) ) ) , sum(S,T), step(C,Z), D = 1..Z, 



T-L+1 #sum [ fail(P) = W 

wait (P,D-1) = W 

true (N) = W 



wlist (S,Q, pos (P) ,W) : not sec (C, pos (P ) ) , 
wlist (S,Q, pos (P) ,W) : sec (Cpos (P) ) , 
wlist (S,Q, neg (N) ,W) ]. 



wait (con junction (S) , D-1 } :- sec (C, pos (con junction (S) )} , set (S, pos (P } ) , 

sec (C,pos (P) ) , wait(P,D-l), step(C,Z), D = 1..Z. 

bot :- sec (Cpos (atom(A) )) , true (atom (A) ) , wait (atom (A) , Z) , step(C,Z). 

% saturate interpretations that are not answer sets 

true (atom (A) } :- elem (atom (A) ) , bot. 
fail (atom (A) ) :- elem (atom (A) ) , bot. 



Listing 3. Disjunctive meta-program (metaD . Ip) for reified extended logic programs. 



bound T-L+1 (or U+1) for weights of literals that do not hold (or hold) is used to indicate 
a violated lower (or upper) bound of the reified #sum constraint. 

Given an interpretation of atoms and the corresponding truth values of further constructs 
in an extended program, the meta-rules in Line 33 and 34 are used to derive bot if the 
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interpretation does not provide us with a supported model. To avoid such a derivation 
of bot, every rule of the reified program must be satisfied, and every true atom must have 
a positive occurrence in the head of some rule whose body holds. 

It remains to check the acyclic derivability of atoms belonging to non-trivial SCCs. To 
this end, the meta-rule in Line 38 determines the number Z of atoms in an SCC labeled C 
as the maximum step at which a fixpoint of T, applied locally to C, is reached. Further- 
more, the meta-rule in Line AO-A 1 derives s c c w ( A ) if the atom referred to by A does not 
have a derivation external to C. (Recall that the positive body elements of rules internally 
connecting an SCC, i.e., rules contributing the SCCs edges to the dependency graph, are 
marked by facts over s c c / 2 ; cf . Listing [T]) The acyclic derivability of atoms indicated by 
sccw (A) is of particular interest in the sequel. In fact, our encoding identifies the com- 
plement of a fixpoint of T in terms of atoms A for which wait ( atom (A) , Z ) is derived. 
To accomplish this, the meta-rule in Line 45 marks all atoms of C as underived at step 0. 
As encoded via the meta-rule in Line 46-47, an atom A stays underived at a later step D if 
there is no external derivation of A (sccw (A) holds) and the bodies B of all component- 
internal supports of A are yet underived at step D-1 (wait (B, D-1 ) holds). The latter is 
checked via the meta-rules in Line 49-52 and 54-55, respectively. The former applies to 
#sum constraints and identifies cases where the weights of literals that do not hold along 
with the ones of yet underived atoms of C exceed T-L, so that the lower bound L is not yet 
established. Similarly, the underivability of a conjunction is recognized via a yet underived 
positive body element internal to the component C. Also note that the falsity of elements 
of C is propagated via the meta-rule in Line 43, so that false atoms, #sum constraints, 
and conjunctions do not contribute to derivations of atoms of C. As mentioned above, the 
complement of a fixpoint of T contains the atoms A such that wait (atom (A) , Z) is 
eventually derived. If any such atom A is true, failure to construct an answer set is indi- 
cated by deriving bot via the meta-rule in Line 57. 

Finally, saturation of interpretations that do not correspond to answer sets is accom- 
plished via the meta-rules in Line 61 and 62 of Listing [3] They make sure that bot is in- 
cluded in an answer set of the meta-program only if it is inevitable wrt every interpretation. 
When considering the encoding part in Listing [3] in isolation, it like meta . Ip describes 
answer sets of a reified program, and bot is derived only if there is no such answer set. 

Our meta-programs meta . Ip and metaD . Ip in Listing |2] and [3] have not yet consid- 
ered facts minimize (J, S) in reified programs, reflecting input ^minimize statements. 
In fact, complex optimization is addressed by the meta-program metaO . Ip, shown in 
Listing |4] It allows for separate optimization criteria per priority level J and weight W 
(in facts wlist (S, Q, E, W) ). Particular criteria can be provided via the user predicate 
optimize (J, W, 0) , where the values card, incl, and pref for refer to minimal- 
ity regarding cardinality, inclusion, and preference JSakama and Inoue 2000l l. respectively, 
among the involved literals E. Such criteria are reflected via instances of cxopt ( J, W, ) , 
derived via the rules in Line 7 and 8-9, where card is taken by default if no criterion is 
provided by the user. At each priority level J, Pareto improvement of a counterexam- 
ple (constructed via the rules in metaD . Ip) over all weights W and criteria such that 
cxopt ( J, W, 0) holds is used for deciding whether a candidate answer set (constructed 
via the rules in meta . Ip) is optimal. To this end, similarity at a priority level J is indicated 
by deriving equal (J) fromequal (J,W,0) over all instances of cxopt (J,W,0) via 
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1 % extract (complex) optimization criteria per priority level and weight 

2 % (relative to user predicate optimize/3; cardinality taken by default; 

3 % Pareto improvement over weights used for comparison at a priority level) 



5 



cxopt (card) . cxopt (incl) . cxopt (pref ) . 



7 cxopt(J,W,0) :- minimize (J, S) , wlist (S,_,_, W) , cxopt (O) , optimize (J, W, 0) . 

8 cxopt (J, W, card) :- minimize (J, S) , wlist (S,_,_, W) , 

9 not optimize (J, W, O) : optimize (J, W, O) . 

11 % verify dominance 

13 equal (J) :- cxopt (J,_,_) , equal(J,W,0) : cxopt ( J, W, 0) . 

15 chain(Jl,J2) :- cxopt (Jl; J2,_,_) , J2 < Jl, 

16 not cxopt (J3,W,0) : cxopt (J3, W, 0) : J2 < J3 : J3 < Jl . 

18 check(J2) :- cxopt (J2 ,_,_) , not chain(Jl,J2) : chain (Jl, J2 ) . 

19 check(J2) :- chain (Jl, J2) , check(Jl), equal(Jl). 



21 bot 

22 bot 

23 bot 



not cxopt(J,W,0) : cxopt ( J, W, O) . 

check ( Jl) , worse ( Jl) . 

check (Jl), equal (Jl), not chain (J1,J2) : chain (J1,J2) 



25 % require non-existence of dominating answer set 

27 : - not bot . 

29 % check inclusion criteria 

31 ndif f (pos (P) ) :- cxopt (J, W, incl) , minimize (J, S) , wlist (S, _, pos (P) , W) , 

32 true (P) . 

33 ndif f (pos (P) ) :- cxopt (J, W, incl) , minimize (J, S) , wlist (S, _, pos (P) , W) , 

34 not hold(P) . 

35 ndif f (neg (N) ) :- cxopt (J, W, incl) , minimize (J, S) , wlist (S, _, neg (N) , W) , 

36 fail (N) . 

37 ndif f (neg (N) ) :- cxopt (J, W, incl) , minimize (J, S) , wlist (S, _, neg (N) , W) , 

38 hold(N) . 

40 equal (J, W, incl) :- cxopt (J, W, incl) , minimize (J, S) , ndiff(E) : wlist (S, _, E, W) . 

42 worse(J) :- cxopt (J, W, incl) , minimize (J, S) , wlist (S, _, pos (P) , W) , 

43 true(P), not hold(P). 

44 worse(J) :- cxopt (J, W, incl) , minimize (J, S) , wlist (S, _, neg (N) , W) , 

45 fail (N) , hold(N) . 

Listing 4. Meta-program for complex optimization (met aO . Ip) on reified logic programs. 



the rule in Line 13. Furthermore, the rules in Line 15-19 are used to chain successive pri- 
ority levels, where a greater level Jl is more significant than its smaller neighbor J2, 
and to signal whether a priority level J2 is taken into account. The latter is the case if 
equal ( Jl ) has been derived at all more significant priority levels Jl. If it turns out that 
a candidate answer set is not refuted by a dominating counterexample, we derive bot via 
the rules in Line 21, 22, and 23: the first rule applies if there are no optimization crite- 
ria at all, the second one checks whether the counterexample is worse (or incomparable), 
as indicated by worse (Jl ) at an inspected priority level Jl, and the third one detects 
lack of Pareto improvement from equality at the lowest priority level. Finally, the integrity 
constraint in Line 27 stipulates bot to hold. Along with saturation (in metaD . Ip), this 
implies that a candidate answer set (constructed via the rules in meta . Ip) is accepted 
only if there is no dominating counterexample, thus selecting exactly the optimal answer 
sets of an input program. The described rules serve the general purpose of identifying un- 
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dominated answer sets, and the remainder of metaO . Ip defines equal ( J, W, 0) and 
worse (J) relative to particular optimization criteria. 

Inclusion-based minimization, indicated via cxopt ( J, W, incl ) , is implemented by 
the rules in Line 31^5. The test for equality, attested by deriving equal (J, W, incl ) 
via the rule in Line 40, is accomplished by checking whether a candidate answer set and 
a (comparable) counterexample agree on all involved literals E; otherwise, ndif f (E) is 
not derived via the rules in Line 31-38. Furthermore, the counterexample is incomparable 
to the candidate answer set if it includes some literal not shared by the latter; in such a case, 
wo r s e ( J ) is derived via the rules in Line 42^3 and 44^5 . In fact, the three C } -minimal 
answer sets of Hi (given in p 1 . Ip), consisting of the rules in (|2]i-(|4]i and the ^minimize 
statement in (|5]i can now be computed in the following way: 

gringo --reify pl.lp | gringo meta.lp metaD.lp metaO.lp \ 
<(echo "optimize ( 1, 1, incl) ." ) - | claspD 

Observe that claspD (' Drescher et al. 2008b . the disjunctive extension of clasp 
JGebser et al. 200Ti) . is used for solving the proper disjunctive ground program obtained 
from gringo. 

In addition to inclusion-based minimization, metaO . Ip implements comparisons wrt 
cardinality dSimons et al. 2002] l and literal preferences JSakama and Inoue 20001) . acti- 
vatable via facts of the form optimize (J, W, card) and optimize (J, W, pref ) 
(along with prefer (El , E2 ) ), respectively. For space reasons, the details are omitted 
here; they can be found in the extended version of this paper ( Gebser et al. 201 ll . The lat- 
ter also provides formal results and arguments demonstrating the correctness of our meta- 
programming technique wrt the specification of optimal answer sets in the background. 

Regarding the computational complexity of tasks that can be addressed using our meta- 
programming approach to optimization, we first note that deciding whether there is an 
optimal answer set is in NP, as the existence of some answer set (decidable by means of 
met a . Ip only) is sufficient for concluding that there is also an optimal one. However, the 
inherent complexity becomes more sensible if we consider the question of whether some 
atom a belongs to an optimal answer set. To decide it, one can augment the reified input 
program (but not the input program itself), met a . Ip, metaD . Ip, and metaO . Ip with 
the integrity constraint : - not hold (atom (a) ) . Then, several complex optimization 
criteria at a single priority level 1 lead to completeness for Sl', the second level of the 
polynomial time hierarchy, thus showing that disjunctive ASP is appropriate to implement 
them. To see this, note that deciding whether an atom a belongs to some answer set of 
a positive disjunctive logic program is El'-complete dEiter and Gottlob 19951 1. When dis- 
junctions fli V • • ■ V flfe in the heads of rules are rewritten to 1 #sum[ai = 1, . . . , a^ = 1], 
the question of whether an atom a belongs to an answer set of the original program can 
be addressed by reifying the rewritten program, adding the integrity constraint : - not 
hold (atom (a) ) , and applying met a . Ip, metaD . Ip, and met aO . Ip wrt several op- 
timization criteria. For one, we can include a ^^minimize statement over all atoms of the 
input program, each associated with a different weight, to exploit the Pareto improvement 
implemented in metaO . Ip for refuting a candidate answer set including a if it does not 
correspond to a minimal model, i.e., an answer set of the original program. Alternatively, 
we can include a ^minimize statement over all atoms of the input program, each hav- 



Complex Optimization in Answer Set Programming 13 

ing the weight 1, and augment the meta-program with the fact optimize (1, 1, incl) . 
We could also use a ^minimize statement over all atoms a,; of the input program along 
with their negation, each having the weight 1, and add the facts optimize (1, l,pref ) 
as well as prefer (neg (atom (ai) ) , pos (atom (a^) ) ) . In view of these reductions, 
we conclude that Pareto efficiency, inclusion, and literal preferences independently capture 
computational tasks located at the second level of the polynomial time hierarchy, and our 
meta-programs allow for addressing them via an extended program along with facts (and 
possibly also integrity constraints) steering optimization relative to its reification. 

5 Applications: A Case Study 

While the approach of Eiter and PoUeres (120061) consists of combining two separate logic 
programs, one for "guessing" and a second one for "checking," into a disjunctive pro- 
gram addressing both tasks, our meta-programming technique applies to a single (reified) 
input program along with complex optimization criteria. In fact, we provide a generic 
implementation of such criteria on top of extended programs encoding solution spaces. 
Hence, our meta-programming technique allows for a convenient representation of rea- 
soning tasks in which testing the optimality of solutions to an underlying problem in NP 
lifts the complexity to Sf -hardness. Respective formalisms include ordinary, parallel, as 
well as prioritized circumscription ( McCarthy 1980 ILifschitz~1 985). minimal consistency- 



based diagnosis dReiter 1987b . and preferred extensions of argumentation frameworks 



(Dung 1995 i. Similarly, Pareto efficiency is an important optimality condition in decision 



making (Chevaleyre et al. 2007 1 and system design (.Gries 2004 J . In the following, we il- 
lustrate the application of our approach on the example of an existing real-world applica- 
tion: repair wrt large gene-regulatory networks dGebser et al. 201()] l. 

Listing|5]shows a simplified version of the repair encoding given in (" Gebser et al. 201()] l. 
It applies to a regulatory network, a directed graph with (partially) labeled edges, repre- 
sented by facts of the predicates vertex/ 1, edge/ 2, and obs_elabel/3, where a la- 
bel S is 1 (activation) or -1 (inhibition). In addition, the data of experiments labeled P are 
provided by facts of the predicates exp/1, inp/2 denoting input vertices (subject to per- 
turbations), and obs_vlabel/3, where a label S is again 1 (increase) or -1 (decrease). 
The regulatory network is consistent with the experiment data if there are total labelings 
of edges and vertices (for each experiment P) such that the label of every non-input ver- 
tex V is explained by the influence of some of its regulators U, where the influence is the 
product S * T of the edge label S and the label T of U (in experiment P). In the practice of 
systems biology, regulatory networks and experiment data often turn out to be mutually in- 
consistent, which makes it highly non-trivial to draw biologically meaningful conclusions 
in an automated way. To address this shortage, several repair operations were devised in 
(IGebser et al. 201()] l. which can be enabled via facts of the form repair (K, J, W) , where 
K indicates a certain kind of admissible repair operations, J a priority level, and W a weight. 
The repair operations R to apply are selected via the rule in Line 14 of Listing |5] and their 
effects are propagated via the rules in Line 18-29, thus obtaining total edge and vertex 
labelings witnessing the reestablishment of consistency. Given that applications of repair 
operations modify a regulatory network or experiment data, we are interested in applying 
few operations only, which is expressed by the ^minimize statement in Line 33. 
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1 % auxiliary concepts 

3 sign(-l;l) . 

4 complement (S, -S) :- sign(S) . 

6 % construct candidate repair 



8 pos (aedge (U, V) , J,w) 

9 pos (eflipCU, V, S) , J,W) 

10 pos (ivert (V) , J,W) 

11 pos (pvert (P, V) , J,W) 

12 pos (vflip (P, V, S) , J, W) 



repair (aedge, J, W) , vertex(U;V), U !- V. 

repair (ef lip, J, W) , obs_elabel (U, V, S) . 

repair {ivert, J, W) , vertex (V) . 

repair (pvert, J, W) , exp(P), vertex (V) . 

repair (vflip, J, w) , obs_vlabel (P, V, S) . 



14 { apply(R) ) :- pos(R,_,_) . 



16 



% construct consistent total labelings 



18 elabel (U, V, S) 

19 elabel (U,V,T) 

20 elabel (U,v, S) 

21 elabel (U,v, S) 

23 vlabel (P,V,S) 

24 vlabel (P,V,T) 

25 vlabel (P,V,S) 



not apply (ef lip (U,v, S) ) , obs_elabel (U, V, S) . 

apply (ef lip (U, V, S) ) , complement (S, T) . 

apply (aedge (U,V) ) , not elabel (U, V, T) , complement (S, T) . 

edge(U,V), not elabel (U, V, T) , complement (S, T) . 

not apply (vflip (P,V, S) ) , obs_vlabel (P, V, S) . 

apply (vflip (P, V, S) ) , complement (S, T) . 

not vlabel (P,V, T) , complement (S, T) , exp(P), vertex (V) 



27 inf(P,V,S*T) :- elabel (U, V, S) , vlabel (P, U, T) , not inp(P,V). 

29 :- vlabel (P,V, S) , not inf(P,V,S), not inp(P,V), not apply (ivert (V) ; pvert (P, V) ) . 

31 % optimize repair 

33 tminimize [ apply (R) = W @ J : pos(R,J,W) ]. 

Listing 5. Encoding of repair wrt regulatory networks and experiment data (repair . Ip). 

A reasonable repair configuration could consist of facts of the following form: 

repair (ivert, Ji, Wi) . admitting to turn vertices into inputs in all experiments. 

repair(eflip,J2,W2) . admitting network modific ations by flipping edge labels . 

repair (pvert, J3, W3) . admitting to turn vertices into inputs in specific experiments. 

repair (vf lip, J4, W4) . admitting data modifications by flipping vertex labels. 

While the kinds of repair referred to by ivert and ef lip operate primarily on a network 
(in view of incompleteness or incorrectness), the ones denoted by pvert and vflip 
mainly address the data (which can be noisy). If we penalize all repair operations uni- 
formly via J = Ji = J2 = J3 = J4 and W = Wi = W2 = W3 = W4, the instantiation of the 
^minimize statement in Line 33 represents ordinary cardinality-based optimization, as- 
sembled in solvers Uke clasp and smodels. However, by adding optimize ( J, W, incl) 
as a fact, we can easily switch to inclusion-based minimization and use a disjunctive solver 
like claspD to solve the more complex problem. While our meta-programs enable such 
a shift of optimization criteria by means of adding just one fact, a direct disjunctive en- 
coding of inclusion-based minimization has been provided in JGebser et al. 201()b : note 
that the latter is by far more involved than the basic repair encoding in Listing |5] Further- 
more, our meta-programming approach allows us to distinguish between different kinds of 
repair operations (without prioritizing them) and optimize wrt Pareto efficiency. To accom- 
plish this, one only needs to pick unequal values for Wi, . . . , W4, where cardinality-based 
minimization wrt each W^ can selectively be replaced by inclusion via providing a fact 
optimize ( J, W^, incl) . Finally, we can choose to rank kinds of repair operations by 
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providing different priority levels Ji, . . . , J4. In this respect, the flexibility gained due 
to meta-programming allows for deploying and comparing different optimization criteria, 
e.g., regarding the accuracy of resulting predictions (cf. (IGebser et al. 2010] l). 

For giving an account of the practical capabilities of our meta-programming approach, 
we empirically compared it to the direct encoding of inclusion-based minimization in 
(IGebser et al. 201 Oi l. To this end, we ran gringo version 3.0.3 and claspD version 1.1 on 
100 instances wrt three kinds of admissible repair operations, resulting in 300 runs each 
with our meta-programs and with the direct encoding. All runs have been performed se- 
quentially on a machine equipped with Intel Xeon E5520 processors and 48 GB main 
memory under Linux, imposing a time limit of 4000 sec per run. To our own surprise, 
more runs were completed in time with the meta-programs than with the direct encoding: 
219 versus ISOQThe disadvantages of the direct encoding show that further gearing would 
be required to improve solving efficiency, which adds to the difficulty of furnishing a func- 
tional saturation-based encoding. In view of this, we conclude that our meta-programming 
approach to complex optimization is an eligible and viable alternative. However, enhance- 
ments of disjunctive ASP solvers boosting its performance would still be desirable. 



6 Discussion 

Our integral approach to modeling complex optimization criteria in ASP brings about a 
number of individual contributions. To begin with, we introduce the reification capaci- 
ties of our grounder gringo along with the associated meta-encoding, paving the way to 
the immediate availability of meta-modeling techniques. In fact, the full version of the 
basic meta-encoding in Listing [T] obtainable at (metasp), covers the complete language 
of gringo, including disjunctions and diverse aggregates. Moreover, our meta-modeling 
techniques provide a general account of saturation and, thus, abolish its compulsory repli- 
cation for expressing complex preferences. Of particular interest is the stability property 
of answer sets serving as implicit counterexamples. Unlike the approach of Liter and 
Polleres (120061 1. our encoding avoids "guessing" level mappings. Also, our target lan- 
guage involves choice rules and #sum constraints (ISimons et al. 20021 . and we are un- 
aware of any pre-existing meta-encoding of corresponding answer sets, neither as candi- 
dates nor as counterexamples. Likewise, related meta-programming approaches for gen- 
erating consequences of logic programs (IFaber and Woltran 2009l l or explanations wrt de- 
bugging queries (lOetsch et al. 201()l l do not consider such aggregates (but disjunctions in 
object programs). 

We exploit the two-dimensionality of ^minimize statements by using levels and 
weights for combining a lexicographic ranking with Pareto efficiency. At each level, 
groups of literals sharing the same weight can be compared wrt inclusion. This is ex- 
tended in (IGebser et al. 201 ll by cardinality-based minimization and the framework of 
(ISakama and Inoue 2000l l. relying on a preference relation among literals (given in ad- 
dition to ^minimize statements); the augmented encoding is also available at (l metasp| l. 
In fact, the approach of Section |4] allows for capturing the special cases of parallel and 



All instances and detailed results are available at jmetasp). 
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prioritized circumscription, investigated by Janhunen and Oikarinen in (120041 I2008I I. An 
interesting future extension is the encoding of optimality conditions for logic programs 
with ordered disjunction jBrewka et al. 2004] l. whose custom-made implementation in the 
prototype psmodels interleaves two smodels oracles for accomplishing a generate-and-test 
approach similar to the idea of our meta-programs. Ultimately, our approach could serve 
as an implementation platform for answer set optimization jBrewka et al. 20031) and the 
preference description language proposed in jBrewka 20041) . Last but not least, our meta- 
programs furnish a rich and readily available source of E^-hard challenge problems, fos- 
tering the development of ASP solvers for disjunctive logic programming. 
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